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Vortex Rings Generated by a Shrouded 
Hartmann-Sprenger Tube 


Jack Wilson 
QSS Group, Inc. 
Cleveland, Ohio 44135 


The pulsed flow emitted from a shrouded Hartmann-Sprenger tube was sampled with high-frequency pressure 
transducers and with laser particle imaging velocimetry, and found to consist of a train of vortices. Thrust and mass 
flow were also monitored using a thrust plate, and orifice, respectively. The tube and shroud lengths were altered to 
give four different operating frequencies. From the data, the radius, velocity, and circulation of the vortex rings was 
obtained. Each frequency corresponded to a different length to diameter ratio of the pulse of air leaving the driver 
shroud. Two of the frequencies had length to diameter ratios below the formation number, and two above. The 
formation number is the value of length to diameter ratio below which the pulse converts to a vortex ring only, and 
above which the pulse becomes a vortex ring plus a trailing jet. A modified version of the slug model of vortex ring 
formation was used to compare the observations with calculated values. Because the flow exit area is an annulus, 
vorticity is shed at both the inner and outer edge of the jet. This results in a reduced circulation compared with the 
value calculated from slug theory accounting only for the outer edge. If the value of circulation obtained from laser 
particle imaging velocimetry is used in the slug model calculation of vortex ring velocity, the agreement is quite 
good. The vortex ring radius, which does not depend on the circulation, agrees well with predictions from the slug 
model. 

Nomenclature 

a radius of vortex ring core 

/ 2 2 

D effective diameter of jet = s - d t 

D s inner diameter of shroud 

d t outer diameter of resonance tube 

E(e) dimensionless energy of vortex ring 
/ frequency 

/(e) dimensionless impulse of vortex ring 
K circulation 

L length of volume (slug) of air in each pulse 

m j et mass flow rate of jet 

impulse mass of air in one pulse 

N formation number 

P vortex ring impulse 

R vortex ring radius 

r distance along radius of vortex ring 

T energy of vortex ring 

t time 

Tjet thrust of jet 

U characteristic jet velocity 

u(t) jet velocity as function of time 
U i average velocity defined in equation (15) 

U 2 average velocity defined in equation (16) 

W vortex ring velocity 

z distance along axis of vortex ring 

c ratio a/R 

f ratio of ring to slug circulation 
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density of air 

duration of one pulse (l/ / ) 
vorticity 
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I. Introduction 

In an experiment to study pulsed ejector thrust augmentation, a shrouded Hartmann- Sprenger tube was used to 
create pulsating flow. 1 A Hartmann- Sprenger tube is a device which, without any moving parts, transforms a steady 
flow into a pulsating flow. By shrouding the device, the flow is channeled into an axially directed pulsed jet, which 
is convenient for ejector studies. This work was performed to establish the properties of this flow. 

It is well-known that an impulsively started flow produces a vortex ring (Elder and deHaas, 2 Das et al. 3 ). Most 
studies of this phenomenon have used a single pulse of water propelled into a tank of water by a piston. The velocity 
and size of the vortex ring can be calculated from the so-called slug model, given the circulation of the flow. The 
circulation is produced at the outer edge of the flow, and is of the order of ULI 2. If the orifice producing the flow is 
an annulus, rather than a circle, one could argue that a quantity of circulation of the order of ULI 2 is produced at the 
outer diameter, and a quantity of equal magnitude, but opposite sign, is produced at the inner diameter. These should 
cancel each other, so that no vortex ring will result. The exit geometry of the shrouded Hartmann- Sprenger tube is 
an annulus, so it is of interest to see whether the train of pulses emerging from it constitutes vortex rings, and if so, 
whether their velocity and size can be calculated from the slug model, which is intended for a single circular pulse. 

In this study, the flow exiting a shrouded Hartmann- Sprenger tube was probed with high-frequency pressure 
transducers, and found to be a train of vortex rings. The probes located the radius of the ring cores, and were also 
used to measure the ring velocities, as well as the velocity of the flow at the exit plane of the Hartmann- Sprenger 
tube. In addition, laser particle imaging velocimetry was used to map the ring velocity and vorticity distribution. The 
measurements are compared with calculations from the slug model. Modifications to the slug model were found to 
be necessary in order to achieve good agreement with the data. 

II. Summary of the Slug Model 

In the slug model, as shown in figure 1, the pulsed jet is assumed to consist of repeated slugs of fluid, each of 
length L , with constant velocity U, and density p, issuing from a tube of diameter D , as would be produced by a 
piston impulsively started with velocity U and stopping after a stroke L. Initially, the model will be in the usual 
form, i.e., for a circular jet orifice, not an annulus. According to Gharib, Rambod, and Shariff, 4 the slugs will 
transform entirely into vortex rings, provided L/D < 4. Above a value of L/D = 4, vortex rings will form with 
circulation, and size, equal to that produced at L/D = 4, followed by a trailing jet containing the rest of the 
circulation. The maximum value of L/D which still produces only a vortex ring with no trailing jet is called the 
formation number, N. Each vortex ring will have a radius R, and a core of “radius” a. Strictly speaking, the core is 
not a circle, therefore does not have a radius. However, it is approximately circular, and what is meant by a , is the 
square root of the core area divided by n. The value of 8 = a/R , is small at small values of L/D , and increases 
monotonically up to a value of about 0.53 at the cut-off point. 

Given the properties of the jet, Linden and Turner 5 have shown that the properties of the vortex, i.e., R , a , and W, 
are uniquely determined as a function of the value of L/D for the jet. To demonstrate this, they balanced the jet 
circulation, impulse and energy with those for the ring, taking the ring values from the work of Norbury. 6 Then with 
values for the jet in the middle of the equations below, and values for the vortex ring on the right hand side, the 
balance equations become 


Circulation: 

K-slug — U L/2 K r i n g — YK s i U g 

(1) 

Impulse: 

P = npUD 2 L/4 = pR 2 K ring /( s) 

(2) 

Energy: 

T = n p U 2 D 2 Ll& = pR K ring 2 E(s) 

(3) 
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8 = a/R 



In these equations, /(c ) and E(z), the dimensionless values of the impulse and energy, are functions of 8 only. 
Following Linden and Turner, these equations can be combined to give 

(L/D)r' -5 = Jtt/2 {p/p R^K ~ g ) / fr/p R K ~ 2 ) (4a) 

= ^fnJJfij/2 IE(e) (4b) 

The right hand side of this equation, and hence L/D also, is a function of e only. Norbury 6 solved the equations for a 
vortex ring numerically, to obtain values of /(c) and E(z ) at a few values of 8. Since the interest here is in small 
values of s, at several values of this ratio, it is preferable to have a continuous (in a) representation. Such a 
representation was found by Fraenkel 7 who used a solution as a series in 8, valid for small 8. Fraenkel’s solutions are 

/(e)=7i(l + 3e 2 /4) (5) 

and 

E(e) = 0.5 In (8/e)- 7/8 + (3 e 2 /l 6) In (8/e) (6) 

However, Fraenkel’ s values do not agree perfectly with those of Norbury, and so Fraenkel’s expressions were 
modified slightly to give better agreement with Norbury’ s values of impulse and energy, by using polynomial fits in 
8, giving 

/(e) =7i (l + 0.303 s + 0.6184 s 2 - 0.3563 s 3 ) (7) 

E(e) = 0.5 In (8/s)- 0.8594 + 0.098 s 2 In (8/s) (8) 

Values of /(c) and E(z) for four values of 8 as given by Norbury, Fraenkel, and the series of equations (7) and (8) are 
given in table 1 . As can be seen, the series of equations (7) and (8) are good fits to the Norbury values, and will be 
used below, and in the plots in figure 2. 

In figure 2(a), L/D is plotted against 8 as a red line, using equations (7) and (8) in equation (4b), with T = 1, as 
assumed by Linden and Turner. Also, results from Linden and Turner are shown as star shaped points, for 8 = 0.2, 
0.4, and 0.6. 

The agreement is quite good, even for 8 = 0.6. The limiting value of L/D = 4 is found at a value of 
8 = 0.55, so higher values of s than this should not be needed. 

The size of the vortex ring is found by combining equations. (1) and (2) leading to an equation for R/D in terms 
of /(s), and hence, through equation (4b), L/D. 
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( 9 ) 


R/D = l/j2n(e)/n =l/ A /2r(l + 0.303 s + 0.6184 s 2 -0.3563 s 3 ) 

Fraenkel also has an equation for the vortex ring velocity, W, namely 

W = (A: n „ g /47ii?|ln(8/s)-l/4-£ 2 {(3/8)ln(8/£)-15/32}] (10) 

For 8 — 0 this agrees with the equation given by Lamb 8 for small core radius. Again this can be modified to agree 
with the calculated results of Norbury better by using the fit 

W = [K ring /47i/?f ln(8/e) - 1/4 - 0. 85 e 2 {(3/8)ln(8/e) - 15/32}] (11) 

Using equations (1) and (9) in equation (1 1), the ring velocity divided by the slug velocity is seen to be a function of 
L/D and 8 only, 

W/U = r 1 5 (L/D)/ 2(1 + 0.303 e + 0.6 1 84s 2 - 0.3563s 3 ) 
x [ 1ii(8/e) - 1/4 - 0.85e 2 {(3/8)ln(8/e) - 1 5/32}] /8tc 


Table 1. Values of the nondimensional vortex ring impulse I(z),e nergy E(z), and velocity W/U 


8 


0.2 

0.4 

0.6 

0.8 

/(s) 

Norbury 

3.226 

3.422 

3.653 

3.888 


Fraenkel 

3.236 

3.519 

3.990 

4.650 


Fit (Eq. 7) 

3.229 

3.419 

3.656 

3.888 

E ( £ ) 

Norbury 

0.9999 

0.6733 

0.5264 

0.4366 


Fraenkel 

0.9971 

0.7127 

0.5950 

0.5526 


Fit (Eq. 8) 

0.9995 

0.6854 

0.5271 

0.4363 

W !U 

Norbury 

0.4363 

0.5353 

0.6019 

0.6522 


Fraenkel 

0.4393 

0.5188 

0.5761 

0.6026 


Fit (Eq. 17) 

0.4381 

0.5272 

0.6034 

0.6516 


and hence can be expressed in terms of L/D. The values of W/U from Norbury, Fraenkel, and equation (12) are also 
given in table 1, again showing good agreement between the Norbury values and the modified fits. R/D from 
equation (9) (green line) and W/U from equation (12) (blue line), again for T = 1, are plotted against L/D in 
figure 2(a). W/U agrees well with results from Linden and Turner (the diamond shaped points). However, the 
calculated values of R/D do not agree with experimental values of R/D from Liess and Didden 9 (the square points) in 
their dependence on L/D. The calculated values decrease with increasing L/D , whereas the experimental values 
increase. 

Shariff and Leonard 10 have plotted ring circulation divided by slug circulation, i.e., T, against L/D , using data 
from Didden, 11 and Maxworthy, 12 and find that this ratio is not unity, and decreases with increasing L/D. They quote 
two expressions for T, but neither fits all the data. Both are close at L/D = 1, with T = 1.4, which does give the 
experimental value of R/D found by Liess and Didden. By finding a value of T that also fits the Liess and Didden 
value of R/D at L/D = 2.5, the constants in the following expression were found, 

T = K nng /K slug = 0.66 + 0.74 / (l/D) (13) 

This agrees with the two expressions quoted by Shariff and Leonard at L/D = 1, and falls between them at L/D = 2.5, 
so it does not seem unreasonable. Admittedly, Shariff and Leonard point out that the reason for the discrepancy 
between the two expressions may be due to a Reynolds Number dependence. This will be ignored here in the 
interests of simplicity. Using T from equation (13), 8, R/D, and W/U were calculated as functions of L/D, and are 
plotted in figure 2(b). Now there is excellent agreement between the calculated values of R/D and the experimental 
values of Liess and Didden. 
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Fig. 2 (a) Calculations of £, R/D and W/U against L/D, using 
the slug model with r = 1 , and, (b) with r = 0.66 + 0.74/(L/D). 
The calculated c is from equation (4b) with equations (7) and 

(8), plotted as the solid line ( ) for (U 2 /U 1 ) 2 = 1, and as a 

dot-dash line ( ) for (U 2 /U 1 ) 2 = 1.1, W/U is from equation 

(12), plotted as a dash line ( ), and R/D is from equation 

12, plotted as a dotted line ( ). Also plotted are calculated 

results from Linden and Turner 5 for s (*), and W/U (♦), and 
experimental values of R/D from Liess and Didden 9 (■). 


It was implicitly assumed above that U was a constant independent of time. This is rarely true in practice. If the 

actual temporal distribution of velocity from the source is represented by u(t) then it follows that L= ( u(t ) dt , and 

Jo 

the jet circulation, impulse and energy are 


Kslug = { X u(tfdt/2 = U ] L/2 

J O 


(14) 
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(15) 


: n p D 2 1/4 ( J \(tf dt J \(t)dt j = n p D 2 L U x /4 


T = n p D 2 1/8 J \(tf dt J 1 u(t)dt = n p D 2 L U 2 2 /8 


(16) 


where the integrals inside the parentheses define Ui and U 2 , and are taken over one pulse. The right hand side of 
equation (4) is now multiplied by (U 2 /U/) 2 . The effect on the curve of e versus L/D of having a time dependant flow 


is shown in figure 2(b), for a case with 

(u 2 /uo 2 = u. 


In summary, given the slug velocity, including 
its temporal distribution, length and diameter, the 
vortex ring circulation, radius, and core radius, and 
hence velocity, can all be determined, but it is 
essential to input T as a function of L/D. Whether 
the circulation for an annular jet can be derived 
from equation (1) will be determined later. 

III. Experimental Apparatus 
A. Pulsed Jet Source 

The Hartmann- Sprenger, or resonance, tube 
consists of a steady, sonic or supersonic jet, which 
is flowing into a closed tube. Under certain 
conditions, a periodic cycle is established in which 
the jet first fills the tube, then a hammer shock 
inside the tube empties the tube, deflecting the jet 
from the tube in the process. When the tube 
pressure has fallen sufficiently, the cycle can begin 
again. What was not known at the start of this 
effort was whether this phenomenon could also 
produce a directed, pulsed jet. To attempt this, a 
cylindrical shroud was placed around the tube and 
jet, to collect the air leaving the tube and direct it 
out the back of the shroud. Hartmann- Sprenger 
tubes have been shrouded previously (Brocher and 
Pinna 13 ), but with acoustic horns, with the 
objective of amplifying the sound. These acoustic 
horns were closed at the end where the source is 
located, and increased in area with distance away 
from the source. The flow from the horn would 
therefore be diverging, and the velocity at the exit 
would be reduced in value from that leaving the 
source. In order to create a more concentrated 
flow, a cylindrical shroud was used in the present 
work. This shrouded tube did produce vortex rings 
traveling along the extended axis of the device. 
The frequency of operation is dependant primarily 
on the length of the closed tube, but also on the 
length of the shroud. 

The shrouded tube used at a frequency of 
550 Hz is shown in figure 3. A Mach 2 
axisymmetric nozzle with a 12.7 mm diameter 
throat was aligned with a resonance tube 152 mm 
in length, internal diameter 16.8 mm, and external 


Pulsed flow out 



i 

t 


Steady flow in 

Fig. 3 Schematic drawing of the shrouded Hartmann-Sprenger tube. 
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diameter, d t = 25.4 mm. This was surrounded by a shroud, of internal diameter D s = 50.8 mm. A needle was aligned 
with the axis of the jet to stimulate oscillations, as demonstrated by Brocher. 14 Use of the needle makes the device 
relatively insensitive to the distance between nozzle exit and the tube entrance. A supply of air at a pressure of 7.8 
atmospheres ensured Mach 2 operation in the nozzle when exhausting directly to the atmosphere. The average mass 
flow was measured upstream of the nozzle, using an orifice. To generate a frequency of 1100 Hz, a plug was 
inserted into the resonance tube, effectively shortening it to 76 mm in length. New tubes and shrouds were built for 
frequencies of 275 and 125 Hz. 

B. Measurements 

A thrust plate, 762 mm in diameter, was used to measure the average thrust produced by the pulsating jet. The 
thrust plate was attached to an Omega load cell model LC601-25, having a range of ±111 Newtons, to provide an 
electrical thrust signal. The signal from the load cell was fed to an Agilent model 34401 A averaging multimeter. All 
runs lasted one minute, during which time the voltmeter stored 180 readings, and then displayed the average value. 
The experimental procedure involved making six tests to read the thrust of the jet. The jet thrust, T jeh defined as the 
average of the six test readings, typically measured 39.4±0.5 Newtons. 

In order to measure the properties of the vortex rings produced, namely velocity, radius and core radius, the rings 
were sampled with a pair of fast response transducers (Endevco model 8530C-50), each built into the nose of 
hemisphere-cylinder body, and separated by a distance of 25.4 mm. The probes were inserted into the flow 76 mm 
downstream of the jet exit. The signals from the probes were displayed on an oscilloscope, using dc coupling. The 
pair was moved radially to generate a plot of pulse pressure versus radius for the vortex rings. This information gave 
an indication of the vortex radius and core radius. Then, by spacing the probes equidistantly from the jet axis, at a 
radius at which the probes gave the minimum pressure signals (the minima being sharper than the maxima), but 
separating the probes axially, a measurement of the vortex ring velocity could be made. This was repeated for each 
frequency. For the 275 Hz frequency only, a laser particle image velocimeter (LPIV), as described by John, Paxson, 
and Wemet, 15 was used to give a velocity and vorticity map of the vortex ring. One of the high frequency pressure 
probes was also placed at the jet exit to obtain the jet exit velocity as a function of time. 

IV. Experimental Results 

A. Vortex Ring Probing 

Signals from the high frequency pressure probes sampling the vortex rings at a distance (usually 76 mm) from 
the nozzle exit, are shown in the oscillogram in figure 4(a), taken at a frequency of 550 Hz, As indicated in the 
sketch in figure 4(b), two probes were used, both aligned with the z axis, and spaced 25.4 mm apart. The upper 
oscillogram is from the probe closer to the axis of the jet, and exhibits a rise in pressure with time, followed by a 
decay, returning to atmospheric pressure at about half the period of the pulses. After the signal returned to 
atmospheric, there appeared to be a high frequency oscillation on it, lasting until the next pulse. As the probe was 
moved radially, away from the jet axis, the value of the peak pressure first increased slightly, then decreased to zero 
at a radial position equal to about 1.2 times the jet radius. At distances greater than this, a pulse was seen which was 
a decrease in pressure, with a minimum pressure well below atmospheric, as shown by the lower oscillogram in 
figure 4(a), which is from the probe further from the axis. The values of the pulse extrema, expressed as pressures, 
are plotted as a function of distance from the jet axis in figure 4(b). For a frequency of 550 Hz, the absolute outer 
edge of the disturbance, where there is no longer any signal, is at a radial position of 3.2 times the jet radius, i.e., 
8 1 mm. These measurements were repeated at each frequency used. 

The positive signals are generated when the probe is acting as a pitot probe, with flow directed towards the 
transducer. The negative signals are partly generated when the probe is in the core of the vortex ring, such that the 
flow is in the opposite direction, and the probe is acting more like a static pressure probe. This will become more 
obvious below. 

As stated above, in order to measure the vortex ring velocity W , the probes were separated axially by a distance 
varying between 50 to 100 mm, depending on the frequency, and were each at a radial distance from the axis such 
that they were sampling the pressure minima. Velocities derived from these signals are given in table 2. Each 
velocity is an average of about 15 readings. Some variation in velocity with probe axial separation, and with position 
of the leading probe relative to the nozzle exit, was observed, but was less than the experimental error. Also listed in 
table 2 are the values of the thrust of the jet T jet alone for each frequency. The thrust was not the same at each 
frequency used, although the mass flow, m jeh was. This implies that the jet exit velocity U was not the same for 
each frequency. Since the vortex ring velocity is proportional to U, it is appropriate for comparison to correct the 
vortex ring velocity to a common thrust of 39.4 Newtons by multiplying the observed value of W by 39.4/7}^. The 
corrected values are also given in table 2. 


NASA/CR— 2005-213576 


7 



If) 

in 

CD 



Atmospheric 

pressure 

Atmospheric 

pressure 


Dper trace zero 
)wer trace zero 


Time 


1 4 
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Fig. 4 (a) Oscillograms of high frequency pressure probe data; (b) plot of the pressure 

pulse extrema versus normalized probe position; (o), data forf = 1100 Hz; ( •), fit 

to data for f = 1100 Hz; (■), data for f = 550 Hz; ( ), fit to data for f = 550 Hz; (□), 

data for f = 275 Hz; ( ), fit to data for f = 275 Hz; (*), data for f = 125 Hz; ( ), 

fit to data for f = 125 Hz. 


Table 2. Measured vortex ring velocity, W, and thrust corrected vortex ring velocity 

f= 1100 Hz f = 550 Hz f = 275 Hz f=125Hz 


W from Pitot Probes (m/s) 

73 

77 

83 

77 

W from LPIV 



79 


Thrust (Newtons) 

40.0 

44.7 

39.3 

33.8 

W corrected to 39.4 Newtons 

72 

68 

83/79 

89 


B. Laser Particle Imaging Velocimetry 

Laser particle image velocimetry was used to measure velocities within the vortex ring, for a frequency of 
275 Hz only. The values of axial velocity along a radius at the centre of the ring are plotted in figure 5. The velocity 
of the vortex ring itself is 79 m/s. This velocity occurs at the vortex radius R. The maxima and minima of velocity 
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are at radii of R - a, and R + a respectively, from which R = 33.3 mm, a = 18 mm, and hence e = 0.54. A sketch of 
the vortex ring is on the right hand side of figure 5. The minimum signal from the pressure probes is at a radius of 
43.2 mm, at which radius the velocity is zero. This radius is close to R + all. For greater radii, the velocity over the 
probes is negative, so they are not then acting as total pressure probes. 

The LPIV data provide a map of vorticity co over the entire field of view. The vorticity along the vortex radius is 
also shown in figure 5. The vorticity is zero at a large distance from the vortex, then after a region of small negative 
vorticity, increases as the distance decreases, peaks, then decreases to zero at a radius of 15.8 mm, and goes to a 
negative peak before returning to zero on the axis. The positive peak of vorticity occurs at the vortex ring radius, 
which is also the radius at which the velocity is equal to the vortex ring velocity. The circulation in the vortex ring 
can be determined from the vorticity by performing the integral 

K r ing = j I cd drdz (17) 

over the whole flow field. The result is K ring =13.7 m 2 /s. The integration was performed numerically from large 
radius to zero. The circulation increased as the radius decreased, reaching a maximum of 23.9 m 2 /s _1 at a radius of 
15.8 mm, then decreasing to the value of 13.7 m 2 /s _1 on the centerline. This is a consequence of the flow coming 
from an annulus, so that vorticity shed from the inner diameter is of opposite sign from that generated at the outer 
diameter, thereby reducing the total circulation. However the negative vorticity shed at the inner diameter does not 
cancel that shed at the outer diameter. 



Fig. 5 Velocities ( — ■ — ) on the r-axis through the vortex ring as a function of normalized radial 
position. The velocities come from LPIV data. V c = K/2na = 130 m/s is the velocity at the edge 

of the core. The vorticity (• ) along the r-axis divided by 100 is also plotted. A sketch of the 

vortex ring is given at the right. 
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Time (sec) 

Fig. 6 Flow velocity at the exit of the Hartmann-Sprenger tube as a function of time, for 
f = 275 Hz. Data from the high frequency pressure probe. 

C. Jet Exit Velocity 

By placing one of the high frequency pressure probes used to sample the vortex ring exactly in the exit plane of 
the jet, a signal is obtained which can be related to the jet exit velocity. This was done at a frequency of 275 Hz 
only. The resulting velocity is plotted in figure 6. The flow starts with a weak shock wave (Mach number = 1.25), 
and, after an initial spike, is approximately constant for a while before falling off. An average flow exit velocity can 
be derived from the thrust, namely 


u - Tjet /fft jet ( 18 ) 

this gives a value of 189 m/sec, which is also shown in figure 6, and is in good agreement with the probe results. 
The probe velocities can be integrated to give a value for ( U 2 /Ui ) 2 of 1.1, the value used in figure 2(b) in deriving 
L/D versus s for a jet with velocity that is not constant with time. 

V. Comparison between the Slug Model and Experimental Results 

The experiments were performed at four different frequencies. Since the mass flow was constant at all 
frequencies, the mass per pulse m p is given by 


nip —J'hjet/f ( 19 ) 

in which /is frequency. But m p is also given by the product of pulse volume and density 

m p = p 7i {p} - d} )z/ 4 (20) 

from which it is seen that L/D is inversely proportional to frequency. In referring to the slug flow model, it is clear 
that the appropriate diameter appearing in L/D is twice the square root of the flow area divided by n , so that for the 

Hartmann-Sprenger tube it is D- ^{d} —d}'). In table 3, the values of L/D for each frequency used are given 

together with values of 8 from equation (4) for the two higher frequencies. For the lower frequencies, since the 
vortex probing showed that the vortex rings were the same size for both frequencies, 8 is the same for both, and 
equal to the value measured from the LPIV data, namely 0.54. Gharib et al. 4 showed that above some value of L/D, 
called the formation number, the circulation of the vortex ring remains constant as L/D increases, with the remaining 
circulation going into a trailing jet. In the case of Gharib et al., the limiting value of L/D was 4, consistent with a jet 
velocity constant in time at 8 = 0.55. If this were the case for the present experiments, the vortex rings for 550, 275, 
and 125 Hz, which all have L/D greater than 4, should be the same size. In fact the 550 Hz vortex rings are smaller 
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Table 3. L/D and 8 for each frequency used 

Frequency (Hz) 1100 550 275 125 

L/D 2.5 5 10 22 

8 0.177 0.386 0.54 0.54 


than those at 275 and 125 Hz. This indicates that the formation number appears at a higher value of L/D than 4. 
From equation (4), with (U 2 /Ui) 2 = 1.1, the maximum value of 8 observed, i.e., 0.54, occurs at an L/D of 6.86. This 
will then be taken to be the formation number N, for these vortex rings. Although the value of N is different from the 
value of Gharib et al., 4 the value of 8 at the formation number appears to be the same, suggesting that it is 8 that has 
a limiting value rather than L/D. 

The probe measurements in figure 4(b) show that the pulse minima are at a value of r/R s =1.7 for the 275 Hz, 
and 125 Hz rings, where r is radius, and R s the shroud radius. The LPIV data in figure 5 indicate that this is the 
radius at which the velocity is zero, and is approximately equal to R + all. Assuming then that the probe minima are 
at a radius of R + al 2, and using the values of s in table 3, the value of R + a can be calculated for each frequency, 
and plotted against L/D. The results are shown in figure 7, together with the LPIV value at 275 Hz, and the 
calculated values from the slug model, assuming a formation number of 6.86. The agreement is very good. 

In figure 8, the observed vortex ring velocities, as corrected for thrust, are also plotted against L/D , together with 
the values calculated from the slug model, again assuming a formation number of 6.86. Two calculated values are 
given, one for a value of K ring =13.7 m 2 /s, corresponding to the overall circulation of the vortex ring from the LPIV 
data, and one for K ring = 20 m 2 /s, corresponding to the value from K ring = T U L 12, i.e., equation (1), with T = 0.768 
at 8 = 0.54. Even though the data show some scatter, it is clear that the velocity calculated with K ring = 20 m 2 /s is 
considerably greater than the observed values. It might be thought that the high value of W calculated from the slug 
theory is because the value of F used in equation (1) is too high. However, T also appears in equation (9), the 
equation for the vortex ring radius, so if T were reduced, the good agreement between calculated and observed ring 
radius would be diminished. In fact the discrepancy is because the circulation calculated using equation (1) is based 
on the circulation created at the outer edge of the shroud only. As shown by the LPIV data, there is vorticity of 
opposite sign shed at the inner tube, which reduces the total circulation. Using the circulation derived from the LPIV 
data, the calculated vortex ring velocities are in much better agreement with the data. 



Fig. 7 (R + a)/D against L/D, with (R + a)/D calculated from the slug 

model ( ), an extension at constant value for L/D greater than 

N ( ), and data from pressure probes (■), and from LPIV 

measurements (□). 
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Fig. 8 Vortex ring velocity plotted against L/D; as calculated from the slug 

model with K r j ng = 13.7 m 2 /s ( ), as calculated from the slug model 

with K r j ng = 20 m 2 /s ( ), together with data from pressure probes (■), 

and data from LPIV (□). 


VI. Comments 

If the flow in the jet as it leaves the shroud is approximately uniform radially, there should be approximately 
equal and opposite circulation generated at the shroud and the inner tube in the jet flow. However, as pointed out by 
Maxworthy, 12 the internal flow in the case of a circular jet is not the only contributor to the circulation, a negative 
circulation appears in the induced flow external to the jet, partially canceling that due to the jet. For the annulus, it is 
not clear what is the boundary condition between jet flow at the inner diameter of the jet and the flow centerline. 
There is probably a trapped region of separated flow, creating a complex flow pattern. It is interesting that the peak 
of the negative vorticity is much closer to the centerline than the outer radius of the inner tube. It is also possible in 
the present flow that some circulation exists in the jet before it leaves the shroud. The present measurement of 
velocity did not have sufficient radial resolution to examine this possibility. Certainly the LPIV results show a 
definite net overall circulation, but the exact value would appear to be a field for further investigation. 


VII. Conclusions 

In conclusion, the flow from the shrouded Hartmann- Sprenger tube has been shown to be a train of vortex rings. 
The vortex ring size and velocity have been measured using probes and LPIV, with good agreement between the two 
techniques. From the leveling-off of both ( R + a) and velocity at high L/D , the formation number was found to be 
6.9, at a value of c = 0.54. The slug model predicts the vortex ring size correctly, provided proper account is taken of 
the dependence of the ratio of ring circulation to slug circulation on L/D , and of the temporal velocity distribution in 
the jet in calculating the relationship between c and L/D. The vortex ring velocity is proportional to the vortex ring 
circulation, which is not given correctly for an annulus by the equation K ring = TUL/2 as this ignores the negative 
circulation generated at the inner radius. However, vortex ring velocities calculated with the circulation derived from 
the LPIV measurements were consistent with the experimentally measured vortex ring velocities. When the time 
dependence of the flow is accounted for, the slug model gives a value of N that is also consistent with the 
experiment, for cut-off at s = 0.54. 
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